Carrier scattering from dynamical magneto-conductivity 
in quasi-neutral epitaxial graphene 
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The energy-dependence of the electronic scattering time is probed by Landau level spectroscopy 
in quasi neutral multilayer epitaxial graphene. From the broadening of overlapping Landau levels we 
find that the scattering rate 1/r increases linearly with energy e. This implies a surprising property 
of the Landau level spectrum in graphene - the number of resolved Landau levels remains constant 
with the applied magnetic field. Insights are given about possible scattering mechanisms and carrier 
mobilities in the graphene system investigated. 
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Understanding carrier scattering in graphene [H is one 
of the fundamental but still unclear issues in the physics 
of electronic properties of this noble material. Initial, 
model prediction for graphene 0, 01 done for short- and 
long-range scatterers independent of carrier density as- 
sumed the relaxation rate 1/r that is directly propor- 
tional to carrier energy e. This is a simple consequence 
of the linear density of states with energy, D{e) ~ |e|, in 
graphene which displays the Dirac band structure with 
a constant Fermi velocity vp. The l/r~|e| relation im- 
plies that the low temperature (Boltzmann) conductivity 
a ~ VpT{Ep)D{Ep) of graphene is insensitive to changes 
of carrier concentration n and/or position of the Fermi 
energy Ep ^ Surprisingly, experiments on 

gated graphene flakes show different behavior, the con- 
ductivity varies linearly with carrier density: <y{n) ~ \n\ 



rule may not apply if 



0, i]. The simple r(e) 
dominant scatterers are of the specific resonant type or 
ripples type [1, 01 , as well as for charge impurities js- 
Nevertheless, the unexpected, strong dependence of a{n) 
on carrier concentration remains puzzling [s, 12| and it 
is not fully clarified whether this behavior is intrinsic to 
graphene or has a well-defined extrinsic origin. 

In this Letter we present a positive test of the 1/r ~ |e| 
rule, which we have performed on quasi-neutral layers of 
epitaxial graphene 13l-ll6|. To probe this graphene sys- 
tem with a fixed Fermi energy we profit of the dynamical 
magneto-conductivity method and follow the energy de- 
pendence of the spectral broadening F of inter Landau 
level (LL) transitions measured in our infrared transmis- 
sion experiments. Notably, the optical conductivity re- 
sponse of graphene involves interband excitations across 
the Dirac point and allows us to to determine the pa- 
rameter F over a broad energy range. We have found 



that in the limit of strongly overlapping LLs, which cor- 
responds to the quasi-classical limit in graphene and in 
which the simple relation F — fi/T{e) is justified, this 
spectral broadening becomes F = a • |e|. This linear de- 
pendence on energy is valid at all magnetic fields and 
for all transition indexes (a = 0.026 in the sample in- 
vestigated). This experimentally determined scattering 
rate l/r(e) ~ |e| dependence rules out specific extrin- 
sic scattering mechanisms, such as resonant scatterers or 
charge impurities, in the epitaxial graphene layers which 
are protected from the environment. A straightforward 
but intriguing consequence of our finding is that the num- 
ber of resolved LLs is constant, independent of the mag- 
netic field. Possible consequences of the 1/r ^ |e| rela- 
tion on the carrier mobility in extrinsically unperturbed 
graphene are discussed. 

The studied sample is a standard multilayer epitaxial 
graphene (MEG) grown on the C-terminated surface of 
silicon carbide (4II-SiC[0001]) with intentional thick- 
ness of 50 layers. A significant part of layers displays an 
electronic band structure as that of an isolated graphene 
monolayer, which results from the characteristic rota- 
tional stacking of these graphene sheets [l3| • This is con- 
firmed by micro-Raman spectra measured on our sam- 
ple, which show a single-component 2D band 18, [l9| . 
with some Bernal-stacked residuals on selected loca- 
tions. The sheets studied in this experiment are quasi- 
neutral [20i] , only several layers close to the interface and 
on the surface of MEG become significantly doped (up 
to 1013 cm-2) 0. 

To measure the infrared transmittance, the sample was 
exposed to the radiation of a giobar, which was analyzed 
by a Fourier transform spectrometer and delivered to the 
sample via light-pipe optics. The transmitted light was 
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detected by a composite bolometer which operated at 
T = 2 K and which was placed directly below the sam- 
ple. Measurements were done in a Faraday configuration, 
using a superconducting magnet. The sample is fully 
opaque in the energy range 85-120 meV and also displays 
a relatively weak transmission around 200 meV. Both ab- 
sorption bands are related to phonon-related absorption 
in the SiC substrate. The spectra presented here are ei- 
ther conventional magneto-transmission spectra normal- 
ized by the zero field transmission T{uj, B)/T{uj,0) or 
differential spectra T{uj, B + AB) /T{uj, B - AB). The 
differential technique is more precise and helps to cor- 
rect for possible magnetic field induced changes in the 
response of the bolometer. 

A representative magneto-transmission spectrum at 
_B = 4 T (normalized to the response at -B = 0) of 
our MEG sample is shown in Fig. [IIa),(b). This spec- 
trum is composed of a series of absorption lines which, 
accordingly to previous studies [3, IS, H^-H^], can be 
easily identified as due to optically active transitions 
between Landau levels L±„ of massless Dirac fermions, 
with the characteristic energies: E±n = ±Vi?-y2e?ii?|n|, 
n = 0,1,2... These are the interband transitions: 
L_(„_|_i)— >L„ and L_„— >L„+i, which are reflected in the 
spectra of our quasi-neutral graphene sheets. From 
the characteristic field dependence of these transitions, 
hu = VF^'iefiB{^ + ^/rlTT), we find vp = (1.025 ± 
0.005) X 10^ m.s~^, in agreement with previous estima- 
tions [l3,[23|- The L_i(o)— >-Lo(i) transition disappeared 
from spectra (due to occupation effect at LL filling fac- 
tor 1/ « 6) at fields of 100 mT which implies density 
n « 10^° cm~^ in the studied layers. 

For a more quantitative analysis of the measured 
magneto-transmission T[uj,B), we first note that the ab- 
sorption of light by our graphene sheets is relatively weak 
and therefore: 



r(w, B) = 1 - C • Re[CT,,(cj, B)]/(eoc) 



(1) 



where the proportionality factor C = 7Voff-i(K^-|-3)/(K^-|- 
1) accounts for the effective number NcS of graphene 
sheets in our sample including the effect of possible par- 
tial coverage of the SiC surface, and for multi-reflections 
of light from the SiC substrate with the refractive index 
K = 2.6 lj|. Re[cra;a;(w, B)] is the real part of the dynam- 
ical longitudinal conductivity of a single graphene sheet 
expressed in Kubo-Greenwood formalism 25]: 



Re[iJxx{uj,B)] 



B 



2hi 



MrAfs-fr)SriEr~Es~nw), 

(2) 

where < < 1 stands for the occupation factor 
of the r"^ LL and Sr{E) for a Lorentzian of width F. 
Mr^s — '^'Fl^\r\,\s\±i are velocity operator matrix ele- 
ments in which 7 = 2 if r = or s = and otherwise 
7 = 1. 
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FIG. 1: (color online) Part (a) and (b): Experimentally ob- 
tained transmission spectrum of multilayer epitaxial graphene 
compared to theoretically expected curves, (c) Broadening F 



as derived from our data for L, 



0(-l) 



(0) and L_i(2) 



-2(1) 



transitions (d) DOS in graphene at i? = 4.0 T as calculated 
for the energy-dependent broadening F = 0.026i5. The broad- 
ening of states around zero energy has been set constant to 
0.026iJi. Several observed absorption lines are schematically 
marked by vertical arrows. 



Equation ([2]) is a simplified version of a more general 
Kubo formula (see e.g. Refs. [2^ and 27) which includes 
initial and final spectral functions. Our simplified for- 
mula can be strictly derived for energy-independent and 
purely imaginary self-energies of the initial and final LLs, 
iVg and iF^, and in the limit of relatively narrow lev- 
els/transitions |25], i.e., Fr,s ^ f>^ and F^^s l-E'r,s|- 
On the other hand, since we deal with interband excita- 
tions only, the condition F^^^ <^ l-E^.s — -E-r+Ls-t-i | does 
not have to be fulfilled and adjacent LLs thus may over- 
lap. With these approximations, the broadening F in 
Eq. [2] corresponds to F = F^ -|- F^, and in a view of re- 
markable electron-hole symmetry of bands in graphene, 
also corresponds to F « 2Fr.s. At a given magnetic 
field, the broadening F can be viewed as a function of 
the resonance energy /frequency F = T^fkUs^r), where 
hus^r = \Er\ + \Es\ « 2|ii^r,s|- In further, we also accept 



the possibility that F — F(fta;) is a continues function of 
the energy /frequency. As a matter of fact, the conductiv- 
ity ^ is practically identical if F(fta;r^s) is fixed for each 
transition or T(huj) is a slowly varying function over the 
interval fiiOr^s ± F(fiWr-s-s)- This latter condition as well 
as Fr.s <C and F^.s ^ \Er.s\ are a posteriori fulfilled 
in our experimental data. 

Since we probe here quasi-neutral graphene layers at 
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low temperatures, we assume that the Fermi energy is 
pinned to the O"^ Landau level and therefore fs = 1 and 
/,. = in the case of our interband (r, s ^ 0), — >Lr 
transitions. As for a pair of degenerate transitions, 
L_i — >-Lo and Lq — J'Li, which involve the O"^ LL, we can, 
without loss of generality, set /i = and /o = /-i = 1 
for any magnetic field, since the sum of intensities of 
L_i — ^-Lq and Lq — >'Li transitions is independent of de 
facto partial population of the Lq level < 2) [28| . 
The phenomenologically introduced broadening parame- 
ter r is the main focus of this work. In general, this pa- 
rameter depends on the magnetic field, on the transition 
index, or, as discussed above on frequency uj. Notably, 
the physical interpretation of T is different in the case of 
well-overlapped or well-resolved LLs. 

In the limit of well separated Landau levels, F = 
F^-I-Fs may correspond to a single-particle quantum level 
broadening Fs of the density of states (DOS) [2^. This 
regime has been theoretically explored in the case of long- 
range and short-range disorder, especially for scattering 
centers insensitive to carrier concentration Q . For short- 
range scatterers, Fg was found to be independent of the 
LL index but to increase as -y/B. In the more complex 
case of long-range scattering, F^ may also vary as Vb 
but the broadening of the O"^ LL is found to be signifi- 
cantly (by a factor v^) larger than that of any other lev- 
els. Charged impurities, if present, lead to a narrowing 
of LLs. The width of the LL is then predicted to oscillate 
with the filling factor and shows an overall 1/B field de- 
pendence, F ^ 1/D{Ef), instead of the ^/B field depen- 
dence 3, 53] ■ As concluded by the authors of Ref. 26, 
this excludes charge impurities as the dominant source 
of scattering in MEG, in which ^/B- ty^p e broadening is 
characteristic of inter-LL transitions f20!| . 

In the limit of densely spaced and therefore strongly 
overlapping LLs, attained in graphene at sufficiently high 
energies - see Fig. [TJd) for illustration - the broaden- 
ing parameter F can be interpreted as the inverse of a 
transport relaxation time, T ^ h/r within a factor of 
2. This interpretation can be justified by the classical 
Drude approach to describe the carrier motion in a mag- 
netic field 2^, sOl , in which the magnetic field effects are 
apparent when the cyclotron frequency {uJc = VpcB/e 
for graphene) is comparable to the carrier relaxation rate 
T~^. The transport relaxation time is the only parame- 
ter which accounts for broadening effects in the dynam- 
ical conductivity response of graphene at zero magnetic 
field [3| and very likely in the magneto-dynamical con- 
ductivity response in the classical limit when the de nsity 
of states is weakly modulated by the magnetic field 31 1 . 
Notably, our assumption F — H/t is favored by the 
electron- hole symmetry in graphene and, therefore, in 
neutral sheets likely identical relaxation rates for elec- 
trons and holes at energies +e and — e, respectively. The 
transport relaxation time becomes comparable with the 
single-particle relaxation time if short-range scattering 



dominates [3^, which is, as shown later on, consistent 
with our data and which thus facilitates our interpreta- 
tion. 

To compare with experiment, we re-consider the (rel- 
ative) magneto-transmission spectrum measured at i? = 
4 T in Fig. [TJ To model this spectrum, we assume that 
axx{^, B = 0) = /Ah. This is also the numerical result 
from Eq. [5] in the limit of strongly overlapping levels. 
Then following Eqs. [I]and[l]we fit T(a;, B)/T(a;, 0) with 
two adjustable parameters/functions C and F(aj). In the 
first and simplest approach, we set F to be constant, inde- 
pendent of the transition energy. As shown in Fig [TJa) , 
we then fail to reproduce the high energy part of the 
spectrum, but can well reproduce the two lowest, best- 
separated transitions. We find C = 1.1 and F = 3.9 meV 
by best fitting the shape of the L_i(o) ~^Lo(i) transition. 

However, most of the spectrum consists of overlapping 
transitions (in fact with the exception of the two lowest- 
energy transitions) that can be very well reproduced by 
assuming F = a • huj. The best values are a = 0.026 
and C = 1.3 from the fit shown in Fig. [Ijb). This latter 
value implies a surprisingly low effective number of layers, 
iVoff ss 2, in the present sample. We related this to a par- 
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FIG. 2: (color online) Differential transmission spectra T{B + 
AB)/T{B — AB). The lowest two spectra have been taken 
for AB = 0.05 T, the other with AB = 0.125 T. For clar- 
ity, the subsequent curves have been shifted by a factor of 
0.045. Experimental curves are accompanied by theoretical 
fits, assuming the broadening parameter F — 0.026fta;. Let us 
note that the condition B/AB > 2n must be fulfilled to pre- 
vent significant overlap of adjacent transitions that involve 
n"^ LL. The dashed line shows the simulated curve for dif- 
ferential transmission at B = 10 T and AB = 0.125 T for 
a constant broadening F = 8 meV, which has been set to 
reproduce the L_i(_2) — J-L2(i) line. 
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tial coverage of the SiC surface by graphene sheets and to 
the partial AB-stacking, e.g. into bilayers recently iden- 
tified as a minor component in MEG 3^ 33|. Perhaps 
another reason could be an increased dielectric screening 
inside MEG [sHj, which can reduce the total absorption 
significantly. Remarkably, the F = 0.026huj rule is uni- 
versal, in the sense that it applies to spectra measured at 
any magnetic field within the range of the present mea- 
surements. This is illustrated in Fig. [2] with simulations 
of differential spectra measured in the region of overlap- 
ping n > 1 transitions, at five distinct magnetic fields. 
The calculated traces, which accurately fit the data, are 
obtained from Eqs.[l]and[2]with a common set of param- 
eters {a = 0.026, C = 1.3). Here we note that despite 
a clear overlap between transitions, each of them is still 
well-defined in a sense that the above discussed assump- 
tion of narrow absorption lines and of slowly varying F 
is indeed well satisfied, F = 0.026fta; <C fiw. 

To sum up, the broadening parameter F seems to be 
constant for transitions between well-separated levels, 
but tends to increases linearly with energy when levels 
start to overlap significantly. This behavior remains un- 
changed up to the highest resolved transitions, i.e., those 
which originate from practically merging LLs. This is 
shown in Fig. [ijd), where the DOS has been plotted 
for the experimentally derived parameters. As discussed 
above, in the regime of strongly overlapping levels, the 
broadening of electronic states is governed dominantly 
by the zero-field relaxation rate and the observed linear 
in energy increase of F indicates a l/r ^ |e| dependence. 
Interestingly, the same conclusion can be drawn using a 
simple semi-classical argument: the number of resolved 
transitions in our spectra is roughly magnetic field in- 
dependent, i.e., the onset of LL quantization remains 
always pinned to a certain level (no^lO, see Fig. [2]). 
The energy of this onset moves to higher energies as 
|e| ^ \/B and the cyclotron frequency corresponding 
to this onset thus follows a linear in energy dependence 
~ VpeB/\e\ ^ \e\. If we employ the classical condition 
for this onset, ujcT ~ 1, we conclude l/r ~ |e|. 

Having estabhshed the H/t = a\e\ dependence, let us 
crudely assume that it is not additionally affected by 
changes in the carrier concentration. The Fermi energy 
independent conductivity of "our graphene" is then a = 
— ■ \ (ct ~ 80e^//i for a — 0.026) and we can speculate 



on a low-temperature carrier mobility, /i 



in such a model system at various charge densities. As 
an order of magnitude, we obtain /i w 10® cm^/(V.s) 
for n = 10^'^ cm~^, which is consistent with our recent 
measurements on quasi- neutral graphene (20| . The mo- 
bility falls down to « 10^ cm^/(V.s) for the carrier 
density as high as ti = 10^^ cm^^, what refiects fairly 
well typical mobilities measured for charged, single-layer 
epitaxial gra phe ne on C-face of SiC or interface layer 
of MEG [2l|, |36| . Our hypothetical graphene with n = 
10^^ cm~^ would have a mobility of /i w 10"^ cm^/(V.s), 



in striking agreement with typical values measured for 
high concentration graphene on Si-face of SiC [s^] . 

The 1/T(e) dependence reported here is expected for 
graphene and points towards conventional scattering 
mechanisms in quasi-neutral epitaxial graphene. Apart 
from charge impurities, excluded already earlier in MEG 
due to \/S-broadening of inter-LL transitions [IHl, we 
can also for certain exclude the dominant role of reso- 
nant scatterers, which produce the mid-gap states and 
therefore modify the predicted 1/r ^ |e| relation. Those 
resonant scatters are by definition extrinsic, and maybe 
important for environment exposed graphene, but are 
naturally expected to be absent in our protected layers. 

To complete our discussion of scattering mechanisms 
in MEG, we focus on the two lowest energy transitions 
of Fig. [TJa),(b). These are the two best separated tran- 
sitions in our spectra, apparently related to nearly non- 
overlapping LLs (|n| = 0, 1 and 2) and therefore, their 
behavior can be tentatively analyzed in reference to the- 
ory in Ref. 2. Figure [ijc) shows that the width of both 
transitions scales roughly as ^/B. Moreover, the width 
of the n ^ transitions is roughly the same if not even 
smaller than the width of the n = 1 transition. Accord- 
ing to the theoretical analysis of the dynamical magneto- 
conductivity in the regime of well-separated Landau lev- 
els, this observation points towards short- rather than 
long-range scattering mechanism (both insensitive to car- 
rier energy). Another possibility is the electron-electron 
scattering which also provides a scattering rate that in- 
creases linearly with energy 38] . 

Concluding, we have investigated the dynamical 
magneto-conductivity of quasi-neutral epitaxial graphene 
layers and studied the carrier scattering in these lay- 
ers. We find that the relaxation rate is linear with the 
carrier energy. This dependence has been initially ex- 
pected for graphene, nevertheless, not clearly observed 
yet. Its interesting implication is that the number of 
resolved Landau levels is fixed, independent of the mag- 
netic field. The dominant scattering mechanism in the 
extrinsically unperturbed layers of epitaxial graphene is 
conventional, very likely due to short-range potentials or 
due to electron-electron interactions. 
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